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Abstract 

We show a right unitary transformation approach based on Susskind-Glogower operators that di- 
agonalizes a generalized Dicke Hamiltonian in the field basis and delivers a tridiagonal Hamiltonian 
in the Dicke basis that is diagonalized by a set of orthogonal polynomials satisfying a three-term 
recurrence relation. Our result is used to deliver a closed form analytic time evolution for the case 
of a Jaynes-Cummings-Kerr model and to study the time evolution of the population inversion, 
reduced field entropy and Husimi's Q-function of the field for ensembles of interacting two-level 
systems under a Jaynes-Cummings-Kerr model. 

PACS numbers: 02.30.Ik,03.65.Fd,42.50.Pq 



I. INTRODUCTION 

The Jaynes-Cummings (JC) model is a fundamental building block in quantum optics; it 
describes the interaction of a qubit with a quantum electromagnetic field under long wave 
and rotating wave approximations. It is exactly solvable pQ and has proven useful to describe 
phenomena as Rabi oscillations [2] and collapse and revivals of the atomic inversion [3], just 
to mention a couple; see [I] for a review on the model. If the number of qubits increases, 
the model, known as Dicke or Tavis-Cummings model, shows many-body phenomena in the 
form of a superradiant phase [5]. Dicke model is also exactly solvable [OHH] and is known to 
show super-fluorescence and amplified spontaneous emission; see pj for a recent review on 
the model. 

In recent years, a general Dicke Hamiltonian including quadratic self-interactions on both 
the field and qubit ensemble was introduced to study the effect of the nonlinearities and their 
relation to the Stark shift, in units of H, 



An exact solution through quantum inverse methods involving Bethe anzats was found jTU] . 

The importance of Dicke model and its generalizations lies on its ability to describe 
more than atoms interacting with the quantized field of a cavity; i.e. lasers. For example, 
it may describe open dynamical cavity-QED systems [TT], ion trap systems [12"] . circuit- 
QED systems [131 E], arid Bose-Einstein condensates interacting with classical or quantized 
electromagnetic fields [ToT - fTT] . 

In this brief report, we present an exact solution, up to the roots of a polynomial, to a 
more general Dicke Hamiltonian than that in Eq.Q by including non-indentical nonlinear 
interactions. Our exact solution allows us to recover the well known time evolution of a single 
two-level system, Jaynes-Cummings dynamics, to verify the validity of our approach. Then, 
we study the time evolution of different ensemble sizes under identical values of quadratic 
self-interactions to illustrate the simplicity of our approach and the results it yields; we focus 
on the dynamics of the population inversion of the two-level system ensemble as well as the 
dynamics of the entropy and Q-function of the field. 
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II. THE MODEL 



Let us consider a system composed by an ensemble of N self-interacting two-level systems, 
hereby referred to as qubits, in the presence of a quantized field and a Kerr medium described 
by the Hamiltonian (in units of K) 

H = H + Hi, 

H = ojfdJa + k (a* a) + x(a 2 + a t2 ) + uo q J z + j^Jl, 
Hi=J=(a + tf)(j + + L), (2) 

, 2 



where the A 2 oc (a + term has not been dropped. The qubits ensemble is described 



J+, J- 



2J 2 , 



Jz, J± 



±J Z , 



a, a 



t] = I. The 



by collective Dicke operators satisfying the su(2) algebra, 
while the annihilation and creation operators for a single mode field satisfy 
rest of parameters are defined as usual, the modified frequency of the field and the two 
level transition frequency are given by Uf = uq — x an d oj q , where is the field frequency, 
the second order nonlinearity that arises from the A 2 term by Xi the Kerr medium by k, 
and the dipole-dipole and ensemble-field couplings by £ and g. Classical dynamics for this 
Hamiltonian and reductions have been considered in the context of two-mode Bose-Einstein 
condensates coupled to radiation fields [TSH2"0] . 

Here, we are going to consider a regime of weak nonlinearities in order to use a squeezed 
states basis for the field, described by the transformation, 

f^e^^), (3) 

u f 

to help us get rid of the second order nonlinearity, x- Then, we can introduce a small 
rotation |21|, 



T 2 = e"f + ^ K n ; , < 1. (4) 

We want to emphasize that this small rotation is valid even in the regime where phase 



transitions appear g c = y/Cof (ujq — £) [El E2] • The small rotation has an effect similiar to 
that of the rotating wave approximation and leads to a Dicke Hamiltonian including a Kerr 
medium and dipole-dipole interactions between the qubit ensemble components, 

H = SJ Z + k (a j a) 2 + 7 J 2 + A (aJ + + a ] j^j , (5) 



where we have moved to a frame defined by the total excitation number N = a) a + J z 
rotating at the frequency of the field and defined the parameters 5 = u q — Uf + 2x 2 /uf, 
7 = x/N an d ^ = %9 ( w / — x) ( w / ~~ ^X 2 ) /VNcjf (uj — 2\ 2 + w a (jf). Note that we have 
taken the self-interaction nonlinearities k and £ a couple orders of magnitude smaller than 
the transition frequencies in order to neglect products of couplings and nonlinearities. Also, 
notice that if Hamiltonian ([5| were the starting point, it is trivial to include an Stark effect 
term, [ihJ z . In the special case of equal self- interactions, k — £, without second order 
nonlinearity, x = 0> the problem can be transformed into the iV-atom maser including Kerr 
nonlinearity and Stark shift [TO] . 

In a previous work we implemented a right unitary transformation to solve a quantum 
Landau-Zener problem for a single two-level system [23] . Implementing an extension of such 
approach for the time-independent model at hand we obtain an evolution operator with the 
form 

U(t) = U A (t)U B (t), U x = e- ti ^\ (6) 
where the auxiliary Hamiltonians 

N/2 N/2 

h a = F &*)m\+ E G(j,h)(\j)(j-i\ + \j-i}(j\), (?) 

j=-N/2 j=-N/2+l 

N/2 N/2-j N/2 -N/2+l+j 

h b = e T, F (^)Pk\m+ E E G{j,h) Pk {\j-i){j\ + \j){j-i\), 

j=-N/2+l k=0 j=-N/2+l k=0 

(8) 



where the ket \j) is a Dicke state and pj is the density matrix for the pure state of the field 
with j photons, are diagonalized in the field basis; i.e., they are given in terms of the photon 
number functions 



G(j,n)=X 





N 




n — 


~2 


+ 3 




' N 








+ 1 







3 (3 ~ 1) 



1/2 r 



N 

n + l + — -j 



1/2 



(9) 
(10) 



4 



III. EXACT SOLUTION 



Here we want to present a different and simpler approach. For this reason, let us define 
the right unitary transformation 



N_ 
2 



ff ] = 1, 



-+k-l 



Ff = l-J2 E f* WO'l. Pk = \k)f f (k\. 



(11) 

(12) 
(13) 



Again, the ket \j) is a Dicke state and pj is the density matrix for the pure state of the field 
with j photons. Then, it is possible to write the generalized Hamiltonian ^ as: 

H = TH sc f\ 



(14) 



where the auxiliar Hamiltonian is given by, 

N_ 
2 

H SC = E fU>n)\j)V\ + 



N 
2 



E 9(j,n)(\j)(j-l\ + \j-l)(j\), 



(15) 



where we have used the notation Hsc to bring forward that this Hamiltonian is semi- 
classical-like because it is only expressed in terms of the number operator, 



9 (j, n) = A 



N ( N 



U(t) 



+ 1 -J'(J-I) 



E 



,1/2 



N 

n + I — - — j 



1/2 



2 V2 

It is possible to express the dynamics of this model as the evolution operator 

1 

m! 



-itH 



-itH 



(16) 
(17) 

(18) 



where powers of the form H m are needed. These powers can be obtained by realizing that 
HscT^THsc — Hg C leads to H m = TH m T^. This yields the evolution operator in the 
reduced form 



(j[t) = fe- ltflsc f ] . 



(19) 



The Hamiltonian Hgc is diagonal in the field basis and is symmetric tridiagonal in the 
Dicke basis; thus it is diagonalizable in the angular momentum basis. The eigenvalues of 
this Hamiltonian can be found by the method of minors and are given by the roots of the 
characteristic polinomial 



Pn+i (v) 



" - f i _ y>™ 



1>\ {i') - If ( -y + l,n j /;, v - i (V) 



(20) 



with 



Po(y) = i, 

Pl (v) = v -f ( —,n 



(21) 
(22) 



Pj iy) 



i> - f [ y + 1 - 3,n 



N 



PMM-rhr + 2 -i-"bi-2W. j> 2 ( 23 ) 



The eigenvectors are calculated from the eigenvalue equation for the Hamiltonian and give 



N 
2 



N_ 

2 



E 4%), E i^i 2 = i. 



2 " — 2 



where the amplitudes answer to the following recurrence relations, 

' N 



f ( y,n ! - " 



" I c { n 



V 2 / 2 

0') _L r. (A £\ Jj) I „ (A I 1 ^ Ji) 



o, 



[/ 0'> - ^ 4" + 9 U, n) + g(j + 1, n) cjf^ = 0, 



/1 4* 



2 



N 



fl,n) 4> =0. 



(24) 
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IV. EXAMPLES 

The time evolution given in the previous section accounts for the full dynamics of the 
system and helps calculating any given quantity of interest. As an example, we will focus 
on the time evolution of the reduced density matrix for the field where the initial state is 

JV 

given by a pure state \tp (0)) = E fc 2 = _iv Cj,k\j)f\ k )> 



N 

Pf (0 = E E c v+i-^ c U-k, k Ui, 3 [P + l + j,tjUl k i [ q + l + -,tJ \p) ff (q\. (28) 



The notation Uij (n, t) = ltHsc ) is used to describe the components of the semi 



ltH SC \ 

classical time evolution operator. This allows us to calculate the mean photon number 
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evolution, 



N 

2 °° / N \ ( N \ 

(*))= P c p+^jA+i~k,kUiAp + l + ^,t\Ul k ip + l + -,t\, (29) 

■ , , JV r>=fl \ / \ / 



and in consequence the population inversion (J z (t)} = (N(t = 0)) — (h(t)}. Another inter- 
esting quantity is the purity of the field, 



P(t) = l-Trpj, (30) 

K 

2 00 / N \ 

Tr p 2 = 2_j C P+ l -3J C 1+°- m ^ mC q+l-k,k C p+o-n,n U ,j ip + l+ "^"'M X 

j,k,l,m,n,o=-fP'1=° 

Uo, m (q + o+j,t)U? )k (q + l+j,t) K n (p + +j,t), (31) 

which is a crude indicator of entanglement in a bipartite system [2U 122] and, often, gives 
equivalent information to that obtained from von Neumann entropy, 

(S f (t)) = -Tr[p f (t)\np f (t)}, (32) 

another measure related to bipartite entanglement [26] . 
For a single qubit, 

H = Kn 2 + -o z + A (a<7+ + a) a J) , (33) 

(34) 



the semi- classical Hamiltonian is given by 



. k (n — l) 2 + f AVn . 

H= \ \ r 2 ,l (35) 
A"\M kti — |, 

and it is possible to give a closed form time evolution operator as 

U(t) = fe- lt " sc f\ (36) 
-itiisc _ e -t K [i+2n(h+i)] j nna to{n)t _ i \P(h)a z + 2Ay^a x ] n(h)t 



[ 2 s2(n) 2 J 

^(n) = <5 + «(l-2n), (38) 

fi(n) = yj\p{n)] 2 + 4n\ 2 (39) 

It is trivial to apply the operator (T) to any given initial state ket (bra) and then apply 
the semi-classical exponential. Figure [T] shows the time evolution of the mean population 
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Figure 1. (Color online) Time evolution of the mean population inversion (a,b), reduced field 
entropy (c,d) and Hussimi's Q-function for the field at times equal to half-minimum (e,g) and 
minimum of entropy (f,h) for a single qubit interacting with a quantized electromagnetic field on 
resonance, 5 = 0, under the Jaynes-Cummings model, left column (a,c,e,f), and under a Jaynes- 
Cummings-Kerr model, right column (b,d,g,h). The initial state for both cases is \ip(0)) = — |) 
with a = 5. 

inversion (first row), entropy of the reduced field (second row) and Hussumi's Q-function of 
the field (third row) for a single qubit as given by a Jaynes-Cummings model (left column) 
and a Jaynes-Cummings-Kerr model (right column). Our results are in accordance with 
those in the literature [31 [27] and we can proceed to sample the dynamics of ensembles. 

For an ensemble of qubits, the task of finding a closed form expression for the time 
evolution becomes cumbersome but it is possible to numerically diagonalize the semi-classical 
Hamiltonian and implement the time evolution of any given initial state. As an example, 
we consider the evolution of ensembles of three, Fig. |2j and twenty five, Fig. [3] qubits. The 
information about the particular initial conditions and parameter values can be found in the 
figures and their captions. At the time, it is not our goal to report and in-depth analysis of 
the dynamics of generalized Dicke models but just to present our diagonalization scheme to 
obtain an exact solution via Susskind-Glogower operators. For this reason we will just briefly 
comment some basic characteristics of the dynamics. By considering an initial state given 
by the separable state consisting of a coherent field and the ensemble in its ground state, 
IV'(O)) = | at)/ 1 — N/2), it is possible to see that the Dicke model presents strong collapse 
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and revivals of the population inversion as long as the mean photon number is larger than 
the number of qubits in the system. A clear collapse of the population inversion is seen in 
any case studied here, up to N ~ a 2 . The strength of the oscillations in the population 
inversion diminishes as the number of qubits in the system gets close to the mean photon 
number of the coherent state but they become ever present at smaller times as we get larger 
ensemble sizes for a fixed value of the coherent state parameter. Meanwhile, the purity 
and entropy of such a Dicke model signals an ever-present entangled state between the field 
and the ensemble as the number of qubits gets close to or equal to the mean number of 
photons; i.e., the plots change from strong, well defined, unmodulated dips in the functions 
to a strongly modulated flat-liner close to the value of a mixed reduced density matrix. The 
Q-function for the reduced field behaves as expected. For a <C N, N + 1 well-defined phase 
blobs appear and evolve half of them clock-wise and the other half counter-clock-wise as 
semi-classicatime goes by. The revivals in the population inversion are associated to the 
overlapping of these phase blobs; a stronger revival corresponding to a better overlapping. 

On the other hand, when an interacting ensemble of qubits is considered under Dicke- 
Kerr dynamics, the collapse and revivals of the population inversion are always weak but 
well defined and periodical. Purity and entropy functions point a return to a quasi-separable 
state on the first revival for the cases analysed with the number of qubits less or equal to 
the mean photon number of the field. The mean value of these functions gradually increases 
with time and some dips appear periodically due to the constructive interference of the 
wavefunction components leading to revivals in the population inversion. Under Dicke-Kerr 
dynamics the phase blobs seem heavily defined by the Kerr process and for a = 5 four phase 
blobs appear and two of them evolve clock-wise while the other two do it counter-clock-wise. 
This process produces an overlap of two and two of the phase blobs leading to a weak local 
minima in the purity/entropy but does not register in the population inversion. It is only 
when the four phase blobs overlap that a pronounced local minima and a revival of the 
population inversion appears. 

V. CONCLUSION 

We have considered a general form of a quantum Rabi Hamiltonian including self- 
interactions in both the atomic ensemble and the field. We present a transformation that 
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Figure 2. (Color online) Time evolution of the mean population inversion (a,b), reduced field purity 
(c,d) and Hussimi's Q-function for the field at times equal to half-minimum (e,g) and minimum 
of entropy (f,h) for a quantized electromagnetic field interacting on resonance, 5 = 0, with three 
qubits under the Dicke model, left column (a,c,e,f), and with three interacting qubits under a 
Dicke-Kerr model, right column (b,d,g,h). The initial state for both cases is | -0 (0) ) = |a)/| — |) 
with a = 5. 



removes the quadratic terms of the annihilation/creation operators and another trans- 
formation that acts as a rotating wave approximation in the limit of weak coupling and 
nonlinearities; these transforms yield a generalized Dicke model. As a side result, we extend 
a previous result based on Susskind-Glogower operators that gives the exact dynamics of 
a Jaynes-Cummings model as the product of two evolution operators. Our main result is 
a different and simpler approach involving Susskind-Glogower operators and right unitary 
transformations. These allows us to represent our generalized Dicke model as a transformed 
semi- classical-like Hamiltonian which is diagonal in the field basis and tridiagonal in the 
Dicke basis; thus the diagonalization of this semi-classical Hamiltonian is known up to the 
roots of its characteristic polynomial. The transformed Hamiltonian gives the time evolution 
of the system and provides access to the dynamics of any quantity of interest. 

We use our result to derive a closed analytical form for the time evolution operator of a 
single qubit interacting with a quantized field in the presence of a Kerr medium, a Jaynes- 
Cummings-Kerr model. Also, we present the time evolution of the population inversion, 
reduced field entropy and Husimi's Q-function of the field for ensembles consisting of three 
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Figure 3. (Color online) Time evolution of the mean population inversion (a,b), reduced field 
entropy (c,d) and Hussimi's Q-function for the field at times equal to half-minimum (e,g) and 
minimum of entropy (f,h) for a quantized electromagnetic field interacting on resonance, 5 = 0, 
with twenty five qubits under the Dicke model, left column (a,c,e,f), and with twenty five interacting 
qubits under a Dicke-Kerr model, right column (b,d,g,h). The initial state for both cases is \ip(0)) = 
\a) f \ — t^) with a = 5. 



and twenty-five interacting two-level systems under a Dicke-Kerr model where the interaction 
and Kerr nonlinearity are equal. 

Our approach can follow the dynamics of hundreds and maybe a few thousands of qubits 
in a simple workstation with efficient programming; this is relevant to the description of 
micromasers and the study of fields interacting with Bose-Einstein condensates to give just 
a couple of examples. 
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